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Abstract 

Quasi-one-dimensional quantum antiferromagnets formed by a (i-dimensional 
hypercubic lattice of weakly coupled spin-1/2 antiferromagnetic Heisenberg 
chains are studied by combining exact results in one-dimension and renor- 
malization group analyses of the interchain correlations. It is shown that 
d-dimensional magnetic long-range order develops at zero-temperature for in- 
finitesimal antiferromagnetic or ferromagnetic interchain couplings. In the 
presence of weak bond alternations, the order-disorder transition occurs at a 
finite interchain coupling. Relevances to the lightly doped quantum antifer- 
romagnets and multi-layer quantum Hall systems are discussed. 
PACS numbers: 75.10.Jm, 75.10.-b, 05.30.-d 
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Low- dimensional quantum antiferromagnets (AFM) exhibit many remarkable properties. 
In strictly one-dimension, transitions into ordered states with broken symmetry is absent. 
For the nearest-neighbor Heisenberg spin-S* chains, the low-lying excitations are gapless spin- 
1/2 quanta (spinons) for half-odd-integer S, whereas a finite energy gap exists for integer S 
||IJ. This profound difference is captured in the effective 0(3) nonlinear cr-model (NLoM) 
description by the value of the topological angle (6 = 2irS) in the 1 + 1-dimensional action 

In two-dimensions, spatially isotropic Heisenberg AFM on an unfrustrated lattice is 
proven rigorously to Neel order in the ground state for S > 1 ||. While no such proof 
exists, it is widely believed that it is the case for S = 1/2 as well. Indeed for d > 2, the 
long-wavelength, low-energy physics governing the interactions between the spin-waves in 
the ordered phase can be described by a d + 1-dimensional NLcrM 0. 

In this paper, we study quasi-one-dimensional quantum AFM, and in particular, the 
disorder-order transition associated with the dimensional crossover. Specifically, we con- 
sider S = 1/2 antiferromagnetic (AF) spin chains with Heisenberg symmetry, arranged in a 
^-dimensional hypercubic lattice, and weakly-coupled by AF or ferromagnetic (FM) inter- 
chain exchange couplings J±. Finite intrachain bond alternations are included to study the 
competition between magnetic order and dimerization. Our strategy is as follows. First, 
interchain coupling is considered at a mean-field level ||,[5] in order to treat the important 
correlations that first develop along the strongly coupled chain-direction. The resulting ef- 
fective one-dimensional theory is transformed into the massive Thirring model whose exact 
Bethe ansatz solution |J is used to obtain the static and dynamical quantities. Then, we 
go beyond the mean-field theory, and show that the order parameter fluctuations can be 
described by an anisotropic d + 1-dimensional NLcrM. Renormalization group (RG) analyses 
are carried out to show that in the ground state, cZ-dimensional magnetic long-range order 
occurs for infinitesimal interchain coupling > 0. In the presence of bond alternation, 
long-range order develops when J_l exceeds a finite critical value. Aside from obvious ap- 
plications to real insulating compounds behaving as weakly coupled AF spin-1/2 chains at 

2 



low-temperatures, we will discuss the implications of our results on the magnetic properties 
of underdoped insulating cuprates and point out the relevance to multi-layer quantum Hall 
structures. 

The starting Hamiltonian of the system is given by 

H = J J £[l + 5(-l) 1 ] S hr ■ S m , r + J ± E s *,r • S,, r+M , (1) 

i,r i> r i/" 

where Sj ir is the spin-1/2 operator at lattice site (i, r) with % and r labeling the sites in the 
chain (z) and transverse to the chain (r M ) directions, fi is summed over the z± = 2(d — 1) 
nearest-neighbors in the transverse directions. The intrachain exchange coupling is AF with 
alternating strengths J(l±<5) > 0, whereas the interchain coupling can be either AF ( J± > 0) 
or FM (J_|_ < 0). We are interested in the case where 8, \ J±\/J <C 1. 

The mean-field decoupling of the interchain term in Eq. ([I]) with respect to AF order in 
the z-direction in spin space leads to an effective Hamiltonian, 

h 1d = j£ [i + 5(-iy] s, • s m -h^i-iys;, (2) 

i i 

plus a constant term Hq = z±N s \J±\m^/2. Here iV s is the number of sites along the chain 
and mo = (— l) l (S^) is the staggered magnetization. Eq. (0) describes a ID AFM in a 
self-consistent staggered magnetic field h = z±\ J±\m . 

Next, we perform a standard Jordan- Wigner transformation Sf = vpjipi — 1/2, = 
ipj exp(i7r S}=i ipjipj)- I n terms of the usual left (L) and right (R) moving fermionic fields 
t/>l and ipR, the resulting theory in the continuum limit is given by 



H[ D = I dz 



-iv(ip[d z ip L - ^ R d z ijj R ) + 2gip[ipl i '4j R ipi 



(3) 



It is well known that the values of v and g obtained in the naive continuum limit are not 
correct in the Heisenberg limit. However, a comparison to the exact excitation spectrum 
of the Hamiltonian in Eq. (0)at/i = 5 = O0 leads to v — irJa/2, where a is the lattice 
constant set to unity hereafter. The terms proportional to h and 5 J in Eq. (H) are easily seen 



in their bosonized forms to be relevant operators of dimension x = 1/2 ||. These competing 
(AF order v.s. dimerization) interactions will induce a mass gap (A) with scaling exponent 
1/(2 - x). Thus, A/v oc (h/v, 5) 2/3 . 

Under a global chiral rotation: ipi — * Gxp(i9/2)ipL, ipn ex P( — iO/tyipR with 6 = 
tan _1 5J/A and Aq = h 2 +5 2 J 2 , terms proportional to h and SJ transform into — A Q (ip\ j ipR+ 
yr'Pl)- The resulting Hamiltonian is then identical to the massive Thirring model with bare 
mass Ao and interaction g, which was solved by Bethe ansatz ||. Following Ref. @], we 
obtain g = 2v = n J and the ground state energy gain per site due to m , 

AE = z ± \J ± \m 2 /2 - J(7/10tt 1 / 3 )(A /J) 4/3 . (4) 

The self-consistent value for the staggered magnetization is obtained by minimizing AE 
with respect to m . We find a critical value for the interchain coupling, 

\Jl\ = (l5nV 3 /Uz x )J6 2 ' 3 , (5) 

which separates an AF phase for \Jj_\ > \Jj_\ where the renormalized mass gap A = 
(Uy/3/5ir)z±\J±\ and 
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from a dimerized phase for | Jj_| < \Jj_\ where A dis = (3\/3/n 2 ^ 3 ) JS 2 ^ 3 and m = 0. 

The mean-field theory predicts a Neel temperature Tjy oc z±J± in the ordered phase. 
While this can be correct when the coordination number z± is large, it obviously contradicts 
the Mermin- Wagner theorem, i.e. AF long-range order should not be possible at any finite 
temperature in d = 2. It is thus necessary to go beyond the mean-field theory and include 
the order parameter fluctuations. To this end, we turn to the dynamical spin correlations in 
the ordered phase. Note that the ordering wavevector Q = (Qj_, it) where Qj_ = (ir, vr, . . .) 
for Jj_ > and Qj_ = (0, 0, . . .) for J± < 0. Since the translation symmetry is broken, 
the uniform and the staggered components of the spins are coupled by umklapp scattering 
with momentum transfer q — ► q + Q. The transverse susceptibility in the random phase 
approximation is therefore given by a 2 x 2 matrix relation, 



X (q,w) = X °(q z ^) [l - l^l/(q±)x°fe,^)] _ , (7) 

where q = (q_i_, q z ), and /(q_i_) = Yin exp(z"qi-/i). Using the equations of motion obtained for 
Hid, and the Lorentz-invariance of H[ D valid at low energies u <C J, it is straightforward to 
show that the components of the ID susceptibility Xum Xus-> an d Xs U are entirely determined 
by X^ssi^z-,^) m the long wavelength limit M. The latter has the following form, 

X°ss(Qz, w) = 2 2 2 + M(u 2 - v\l). (8) 

A^ + v z q 2 — uj 2 

Here the pole in the first term arises from the lowest energy triplet excitation which cor- 
responds to an added fermion in the Thirring model. The function A4(x) contains the 
contributions from the continuum involving particle-hole excitations in the Thirring model. 
The latter has a threshold singularity at and a vanishing spectral weight below to = 2 A. 
Thus x(q, uj) in Eq. (0) is dominated by the collective excitations for uj < 2A. We will 
neglect the contributions from A4(x) in Eq. which is equivalent to the single mode 
approximation (SMA). The constant w in Eq. (H) is then fixed at w = A 2 /z±\J±\ by the 

condition X° s (0)0) = 1/ Z ±\J±.\- 

Solving for x(q, uj) in Eq. (0) using the SMA, we obtain the staggered transverse suscep- 
tibility 

, , A 2 l-h 2 /A 2 
«iK±| oj^-uj 2 

where u; q is the gapless spin wave dispersion (Goldstone modes), uj 2 = A 2 (l — f(qj_)/z±) + 



v q 2 in terms of the weakly modified mass gap A = Awl — h 2 /A 2 and velocity v 



vJl — h? / A 2 due to the interchain correlations ||. The uniform static susceptibility, 

h 2 /A 2 1 , . 

^ = z ;\jAi-h 2 ,A 2 - < io > 

In the limit |JJ — > and 5 = 0, x± — 1-07(1/tt 2 J). The close agreement of the latter with 
the exact ID result (l/7r 2 J) suggests that the SMA is a rather accurate description. When 
<5 7^ 0, x± an d thus the spin stiffness vanishes linearly as |Jj_| — > \J±\. 
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The interactions between the AF spin waves can be described by the 0(3) quantum 
NLaM ||. From the susceptibilities derived above, the d + 1-dimensional Euclidean action 
is given by (setting % — 1), 



Here n T (z,r) is a three-component unit-modulus vector field. It represents the local ori- 
entation of the AF order parameter. The discrete sum runs over the neighboring lat- 
tice sites in the transverse directions. As usual, (3 = 1/ksT, and A is a spatial cutoff 
at which the coupling constant g = v/p° s and p° = x±v 2 . The anisotropy is contained 
in R = A 2 /z±v 2 ~ 0.97z±(\J±\/ J) 2 <C 1. In terms of the unit vector field, the trans- 
verse spin susceptibility x s (q, uS) = ml(n + (q, u;)n~(q, u)) and the staggered magnetization 
m = m (n z )\ T=0 . 

The RG analysis of «So is subtle. Let us consider the case when 5 = 0, i.e. for vanishing 
bond alternation. Notice that we did not keep track of the topological term explicitly since 
at 9 = 7r it does not renormalize under the RG. The effect of the latter is however crucial 
for the renormalization of the coupling constant g in the 1 + 1-dimensional sector in the 
limit R — * |TIJ. In the presence of the topological term, g flows to a finite fixed point 
value g(oo) and the correlation length is infinite, whereas g — > oo and the system develops a 
finite correlation length £o-/A ~ e 27T ^ 9 in its absence. Thus, during the RG transformation of 
Eq. flnf) , if the anisotropy is large enough such that V^R/A <^ 1/£ CT , further renormalization 
using the d + 1-dimensional RG cannot eliminate this finite correlation length, the correct 
treatment of the S = 1/2 system must include the effect of the topological term. On the 
other hand, in the opposite limit where 1 ^> V^R/A ^> 1/^, the long-wavelength physics 
is essentially controlled by the d + 1-dimensional RG and the topological term would not 
make a qualitative difference in the ordered phase. Below, we consider the two situations 
separately. 

For \^R <^ e~ 2n ^ 9 , we follow the analysis of dimensional crossover [ [TTJ| , |TT| . Since R is 
exponentially small, the RG in the 1 + 1-dimensional sector can be performed independently 




(11) 
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by integrating out high-momentum modes until the effective couplings become comparable 
in all directions at a larger cutoff Aj_, i.e. when R/A 2 g w l/A\g(Aj_). The <i+ 1-dimensional 
RG is switched on thereafter. For large Aj_, the coupling g(A±) flows towards its limiting 
fixed point value which is of order one. Thus, the crossover length scale is Aj_ ~ A/\/R. 
Since the scaling dimension of the n-field is zero in the ordered phase, there is no need to 
rescale the latter in the interchain term in Eq. (|lT|). Taking the continuum limit in the 
transverse directions by absorbing the cutoff Aj 2 into defining the derivatives, Eq. (|ll|) is 
reduced to a continuous, isotropic action at an isotropic cutoff Aj_, 

Si = ~ n / d d ~'r / dz / dr[(d r n) 2 + (d z n) 2 + («9 T n) 2 ], (12) 

2g d+1 J J Jo 

where g^+i — g(A±)A d ^ 1 is the bare coupling constant for the d + 1-dimensional RG. For 
R — > 0, g{Aj_) — > g{oo). Thus, the stability of the AF ordered state at infinitesimal R is 
determined by whether g^+i is smaller than the critical coupling g d+1 of the d+ 1-dimensional 
RG. The latter predicts a T = fixed point at g c d+1 / A^ 1 = (d - l)2 d 7r d ^T(d/2) to one-loop 
order 

Since the exact solution of Hi D in Eq. (g) describes the fixed point physics in the 1 + 1- 
dimensional sector, g(oo) = \jx\_v can be calculated using the results in Eqs. (Q) and (P~D[) . 
(i) For 5 = 0, g(oo) = b[l - {2/bir)z ± \ J±|/J] 1/2 with b = 810/14vr 2 , which approaches the 
value g(oo) = 5.86 from below as |Jj_| — > 0. Thus, the condition for AF order, g~d+i = 
gd+i(Aj_)/g^ +1 < 1, is satisfied and improved further with increasing £j_|J|_|. This one-loop 
result is consistent with the numerical series expansion analysis in d = 2 for AF coupled 
chains [|1(|. We therefore conclude that in the absence of bond alternation, long-range 



order develops for infinitesimal AF or FM interchain couplings in d > 2. The physical 
origin of this behavior should be traced back to the gapless power-law correlations in the 
spin-1/2 Heisenberg chain. For d = 2, a finite temperature fixed point does not exist in 
the NLcrM, the ordered phase is stable only at T = 0. This corrects the naive mean-field 
prediction of a finite T N . (ii) When 6 is finite, \J]_\ ^ in Eq. (^|). For z±\J±\/J <^ 1, 
we find g(oc) ~ 6(1 — \ J±/ J±\ 3 )" 1 which diverges as |Jj_| is reduced toward \J±\, indicating 



a transition into a disordered phase with finite dimerization. The critical coupling is at 
a \Jj_\ (> \Jj_\ predicted in the mean-field theory) where g~d+i = 1. To one-loop order, 
\Jj_\ = \J±\(1 — #J+i) -1 ^ 3 - The renormalized spin stiffness vanishes on approaching the 
transition according to p s = p°(l — (jd+i)- Since {{n z ) 2 ) = (1 — (jd+i), the renormalized 
staggered magnetization m 2 = J71q(1 — g~d+i)- 

In the above discussion, the coupling constant g defined at cutoff A in Eq. (|TT|) was taken 
to be close to the (order one) fixed point value of H±d- From the point of view of the effective 
NLcrM, this does not have to be the case. If g <C 1, the independent 1 + 1-dimensional RG 
should be replaced by the d + 1-dimensional RG once the reduced anisotropy is in the range 



1 > VR' > e~ 27r/9(A,) at a cutoff A' = A^jR'g/Rg(A'). Further analysis then belongs to the 
second parameter regime which is also relevant for a large-S AFM where g ~ 2/S. 

For 1 3> VR ^> e~ 27r / 9 , because R is no longer exponentially small, the 1 + 1-dimensional 
RG does not transform independently. Instead, one should treat the d + 1-dimensional RG 
with anisotropic cutoffs and expect the residual anisotropy to persist down to the d + 1- 
dimensional fixed point. Taking the continuum limit in the transverse directions at cutoff 
Aj_ = A/\^R in Eq. ([11]) leads to an action like the one in Eq. (|12|) with bare coupling 
gd+i = gA d ^ 1 . The important difference is that the momentum cutoffs are now anisotropic, 
|kj_| < k™ = 7r/Aj_, \k z \ < k™ = 7r/A. We have carried out the momentum shell RG by 
integrating out modes in the high-momentum layers of the ci-dimensional box in fc-space 
with k™ > |kj_| > k™e~ l and k™ > \k z \ > k™e~ l , where e l is the length rescaling factor. At 
T = 0, the one-loop RG equation is dgd+i/dl = (1 — d)gd+i + g^ +l / g d+ i{R)- The critical 
coupling g d+ i(R) depends on the anisotropy R. For d = 2, we find, 



g c 2+x {R) = 2ttA ± VR 



ln(Vl + R + v^R) + y/R\n(Jl + 1/R + 



-i 



(13) 



For R <C 1, g^+iiR) ~ 27rA_ L /ln(2/v A R)- The ratio of the bare coupling to the fixed 
point value is thus §2+1 (R) ~ (g / 2tt) \n(2 / \/R) <C 1 in this parameter regime. In fact, for 
all d > 2 one can show g d+1 (R) oc tcA 1 ^ 1 / ln(2/\AR) such that ga+i(R) 1. Thus the 
conclusion is again that the weakly-coupled spin chains are in the ordered phase. When 



(5^0, the uniform susceptibility x± decreases with decreasing |Jj_|. Thus g increases 
until y/~R ~ e~ 27V / g . The ordered phase becomes unstable at a critical value R c defined by 
gd+i(R c ) = 1 where a transition into the dimerized phase takes place. Close to the transition, 
p s = p° s y/R{l - 9d+i(R)) and m 2 = m§(l - g d+1 (R)). 

The finite temperature properties close to the transition are described by the scaling 
behavior of the d + 1-dimensional quantum NLoM. The correlation length (£) at low- 
temperatures can be obtained by integrating the one-loop RG equations ||. For d = 2, in the 
quantum critical regime, £ = {hv / 7cbT)27tA/ g 2 +i(-R c ) at #2+1 (-Rc) = 1, where v is the effec- 
tive spin-wave velocity along the chain direction. For R c <C 1, £ oc (/u)//cbT)\/RQ hiy/R^\. 
In the renormalized classical regime, £ oc \fR\ In \fR\{Tiv /ksT) exp (2irp s /kBT). Thus the 
effects of anisotropy enter the prefactors of these universal functions. 

Recently, it was proposed that the destruction of AF long-range order in the lightly- 
doped insulating cuprates may be explained by the dimensional crossover in an effective 
low-energy theory of a 2D Heisenberg AFM with increasing anisotropy ||12j| . The important 
question is whether the effective spin quantum number of the AFM corresponds to an integer 
or a half-odd-integer. Since the effective theory is motivated by the physics of the striped 
phases [P^-[To| with a large correlation length compared the distance between the stripes, the 
effective spin presumably depends on the number of spin-1/2 chains between the neighboring 
stripes. The latter is, unfortunately, not known precisely in the lightly-doped regime. It 



was shown in Ref. ||12|| , that for integer effective spins (appropriate for the case of an even 
number of spins between the stripes), AF long-range order disappears below a finite critical 
interchain coupling related to Eq. (|T3|). In the case where the effective spin corresponds 



to a half-integer, our results show that AF long-range order is stable for arbitrarily small 



interchain couplings in agreement with previous theoretical work In order to destroy 
AF order in this case, small dimerization is necessary in addition to anisotropy. Note that 
the amount of dimerization (bond alternation) necessary to disorder the system and open 
up a spin gap is S c oc (J±/J) 3 ^ 2 which can be very small, whereas S c ~ 0.3 on a dimerized 



square lattice Heisenberg AFM without spatial anisotropy |jl6 



The properties of the S = 1/2 AF spin chains can provide useful insights into the integer 
quantum Hall transitions. In a single layer quantum Hall structure (QHS), the latter is in 
the universality class of the dimerization (spin-Peierls) transition of an SU(2n) AF quantum 
spin chain in the limit n — > [[H]]. The double-layer (or spin Landau level mixing) case 
corresponds to two FM coupled spin chains These results suggest qualitatively similar 
phase structures in the n = 1 and n — > cases. The difference in the universal properties 
of the phase transitions can be summarized by the changes in the critical exponents with 
n. A large number of coupled multi-layer QHS corresponds naturally to an array of FM 
coupled spin chains, where the interchain coupling originates from tunneling between the 
layers. The quasi-lD spin chains may order for infinitesimal |Jj_| > 0, as we have shown 
for n = 1. The resulting spin wave spectrum just describes the diffusive modes that would 
appear in a disordered metal, suggesting the formation of a metallic phase between the 
insulator/quantum Hall states, consistent with recent numerical simulations fl9|| . The cor- 
responding phase transitions are thus in the universality class of the U(2n)/U(n) x U(n)\ n ^o 
NLcrM. Interestingly, the latter also describes the 3D Anderson transition in the presence of 
time-reversal symmetry breaking. 

The author thanks A. H. Castro Neto for a discussion. This work was supported in part 
by an award from Research Corporation. 
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